Study of Percolative Transitions with First-Order Characteristics 
in the Context of CMR Manganites 
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The unusual magneto-transport properties of manganites are widely believed to be caused by 
mixed-phase tendencies and concomitant percolative processes. However, dramatic deviations from 
"standard" percolation have been unveiled experimentally. Here, a semi-phenomenological descrip- 
tion of Mn oxides is proposed based on coexisting clusters with smooth surfaces, as suggested by 
Monte Carlo simulations of realistic models for manganites, also briefly discussed here. The present 
approach produces fairly abrupt percolative transitions and even first-order discontinuities, in agree- 
ment with experiments. These transitions may describe the percolation that occurs after magnetic 
fields align the randomly oriented ferromagnetic clusters believed to exist above the Curie temper- 
ature in Mn oxides. In this respect, part of the manganite phenomenology could belong to a new 
class of percolative processes triggered by phase competition and correlations. 

PACS numbers: 71.10.-w, 75.10.-b, 75.30.Kz 
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I. INTRODUCTION 

The understanding of the self-organized behavior 
of transition-metal oxides, such as manganites and 
cuprates, has developed into one of the dominant sci- 
entific themes of condensed matter physics. Phase co- 
existence and competition between different kinds of or- 
ders involving charge, orbital, lattice, and spin degrees 
of freedom leads to physical "complexity" as a charac- 
teristic of their behavior. The balance between compet- 
ing phases is subtle, and small changes in composition, 
magnetic fields, or temperatures lead to large changes in 
the material properties. In the context of colossal mag- 
netoresistance (CMR), these phenomena motivatjed the 
nanoscale phase separation theory of manganitead, and 
its concomitant percolative transition from the metal to 
the insulator. Many experiments have reported results 
compaJjihle with phase separation and percolation in Mn 
oxidesHElo. Recently, CMR has been observed in simula- 
tions of simple models of competing phases, between Tq, 
the ordering temperature, and T*>Tq, the temperatute 
where coexisting clusters start developing upon coolingB. 
In the presence of magnetic fields the preformed ferro- 
magnetic (FM) clusters align their magnetic moments, 
and they percolate. This encouraging agreement between 
theory and experiments suggests that the basic concept 
behind the CMR phenomenon -phase competition- has 
been unveiled. 

However, several experimental facts still await for a ra- 
tionalization within the phase separation scenario. Noto- 
rious among them is the abrupt character of the metal- 
insulator transition (MIT) in some manganites, result in- 
compatible with standard percolative phenomena where 
the transition is smooth. The results in Fig.l illustrate 
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FIG. 1: Results illustrating the "abrupt" character of the 
metal-insulator transition in some manganites (compositions 
indicated), (a) Dependence of residual pdc on chemical sub- 
stitution at (•) and 4 kOe (O)0. The solid line represents 
the prediction of standard percolation, while the dotted one is 
a guide to the eye. The A are results of a simulation using a 
broad distribution of bonds (q=0.857|^. (b) Resistivity vs. T, 
at magnetic fields (•) and 7 T (Opp. (c) Field dependence 
of the spin-wave stiffness at r=40 Ku. 



with y, with a conductivity exponept t much larger than 
predicted by standard percolationnl. (b) Far more chal- 
lenging, Pdc of Smo.55Sro.45Mn03 changes from insulator 
to metal, in a first-order transitionQ. Although most Mn 
oxides have continuous MITs, their pdc vs. T curves of- 
ten show rapid changes and hysteresis loops, in the same 
regime where cluster percolation seems to occur, (c) 
Varying magnetic fieids, Pro.yCao.sMnOa also presents 
a first-order-like MITlZl. These abrupt transitions cannot 
be rationalized within standard percolation theorycl, but 
they are foiind in regimes that exhibit some aspects of 
percolationcl. 

The theory recently presented by Burgy et alJi relies 
on a CMR state above Tc made out of preformed FM 
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clusters, with random magnetic moment orientation (a 
similar picture was also qualitatively discussed by Ue- 
hara et alfl). This inhomogeneous state is stable due to 
the presence of quenched disorder. The random orienta- 
tion of the cluster magnetic moments can be achieved by 
having "walls" between them made out of the compet- 
ing insulating phase. Otherwise, when two FM clusters 
with different orientations become in contact they will 
align their moments. An insulating wall can prevent this 
process, by drastically reducing the interaction between 
domains. However, when a relatively small magnetic field 
is switched on, these preformed large FM moments can 
easily rotate into the direction of that field. When this 
occurs the insulating walls are no longer needed to pre- 
vent moment alignment, and they melt. This last process 
leads to a percolation, namely, as the FM clusters rotate 
their moments they also slightly increase their size to oc- 
cupy the space left by thejnsulating walls, allowing for 
charge conduction to occum. Although this last step ap- 
pears percolative it is likely non-standard and abrupt, 
since it involves large rounded clusters that are not very 
ramified (contrary to standard percolation), due to sur- 
face tension effects which are of relevance in a regime of 
phase competition. 

It is the main purpose of the present paper to describe, 
phenomenologically, this percolation of smooth objects, 
believed to describe CMR manganites. It will be argued 
that first-order characteristics can be found in percola- 
tive processes where surface tension effects are relevant. 
It is important to clarify that the full description of the 
CMR process is far more involved, with the random dis- 
tribution and subsequent field-induced alignment of FM 
moments playing a key rolcir. The first steps in that direc- 
tion were already describedB. Here our goal is the follow- 
ing: if the FM clusters are assumed already aligned (and 
in that case the spin degree of freedom can be dropped) , 
how does the subsequent percolation occurs? It will be 
argued that this transition can be fairly abrupt, even 
discontinuous. The paper is organized as follows. In Sec- 
tion II a variety of non-standard forms of percolation are 
described, searching for situations where first-order tran- 
sitions are induced. The main conclusion of this section 
will be that achieving a discontinuous behavior is highly 
nontrivial and difficult. In Section III, a novel solution is 
proposed involving deterministic rules within a cellular 
automaton system. A discussion of results is presented 
in Section IV, while Section V contains the main conclu- 
sions. 



II. FORMS OF PERCOLATION 

Standard percolation (SP) is a successful statistical 
model for a variety of heterogeneous systemsll. It has 
been used to describe diffusion through porous media, 
conduction through random networks, and forest fires, 
among others. In spite of its successes, however, SP has 
several shortcomings, including the failure to account for 



correlation effects. Indeed, in SP each site (or bond) is 
occupied independently of its neighbors. This feature is 
responsible for the ramified structure of the usual per- 
colation clusters. In Mn grides the coexistence of large 
clusters has been reportedaQ, but the cluster distribution 
-while complex and involving many lengths scales- is not 
of the very ramified form characteristic of SP. A possible 
explanation lies in the long-range character jpf the elas- 
tic distortions induced by holes in Mn oxidesO, as well as 
surface tension effects as described in the next section. 
These factors provide a natural tendency toward smooth 
surfaces of coexisting metal-insulator clusters. In addi- 
tion, in physical systems with phase separation, there is a 
local tendency to homogenize the landscape, minimizing 
the interface tension. 

While simulations of models considering the elastic ef- 
fects through cooperative Jahn- Teller phonons are CPU 
time expensive, the cluster rounding effects appear even, 
in simpler Hamiltonians, such as the one-orbital modeltl 
(the actual form of this Hamiltonian has been extensively 
discussed in previous literature, and it will not be repro- 
duced here). Monte Carlo (MC) results in Fig.2(a,b) il- 
lustrate this phenomenon. In the simulations, the Heisen- 
berg coupling Jaf between localized t2g spins is selected 
such that some regions of the studied lattice are in a 
FM state if they were isolated, while others present some 
cormpeting spin arrangement (antiferromagnetic (AF) or 
fiuxcl). However, when the FM and non-FM regions are 
in contact, there is a considerable influence from one an- 
other, and the system tends to homogenize its trans- 
port properties. For example, the ladder in Fig. 2(a) 
would prefer to be in a perfectly AF spin state if iso- 
lated. However, by immersing the ladder into a ferro- 
magnet, the spins cant appreciably, allowing for substan- 
tial charge hopping across the ladder. A thin AF re- 
gion cannot survive the "pressure" of a FM environment, 
and it transforms into a fairly good conductor. Fig. 2(b) 
shows similar evidence, now using competing FM and 
flux states. A finger-like region emerging from the lower- 
half, with coupling favoring the flux phase, also does not 
survive in its original form but spins rotate to allow for 
conductiont2l. Modeling percolative manganites with the 
metal/insulating character of a link independent of its 
neighbors can only be a crude approximation, and it may 
cause the present discrepancy between theory and exper- 
iment. If the clusters could be made more 'rounded' the 
transition would be more abrupt, and hopefully in better 
agreement with real manganites. The exponent t may 
depend on the short-scale real-space metal-insulator ar- 
rangement, and changing the typical cluster shapes may 
alter t as well. 

The physics of manganese oxides, in particular the 
abrupt change in resistivity as a function of both field 
and temperature, is supposed to emerge from the re- 
alistic iQne- and two-orbital models that are widely 
studiedcl (supplemented by the addition of quenched dis- 
order). However, simulations using realistic Hamilto- 
nians for manganites -containing the rounding effect 
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FIG. 2: (a) Typical snapshot of an 8x8 MC simulation of 
the one-orbital model (Hund coupling Jh=cx3, T=0.01t, t the 
hopping amplitude, and density 0.8). Jaf is everywhere, 
but in the bonds of the ladder where they take the value 
0.12t. The average spin-spin correlation in the ladder is -0.4 
(effective "double-exchange" hopping 0.54), as opposed to - 
1.0 (effective hopping 0) when isolated from the FM region 
(effective hopping 1). (b) Similar to (a) but for density 0.5. 
Here Jaf=0 is in the upper portion minus the firiser-like 
region. In the rest JAF=0.15t, favoring a flux stateO. The 
effective hopping in the flux region is 0.7, while in the finger it 
is as higlf-^s 0.9. (c) The different neighborhoods in Vichniac's 
notationt3. (d) Typical configuration for standard percolation 
on a 128^ lattice at (or just below) the critical density, (e) 
Same as (d) but for correlated disorder (7=0.2). (f) Same 
as (d) but for the annealing rule (MGE5). (g) Same as (d) 
but for the nucleating rule (MGE4). The latter is particularly 
promising since a small increase in p leads to a fully saturated 
final state, in a first-order transition. 



spontaneously- are simply prohibitive in the percolative 
regime if critical exponents are needed, and one must rely 
on semi-phenomenological approaches, as in the present 
paper. Even if we accept this simplification, there is, a 
priori, no obvious way to introduce the proposed local 
cluster compactification in percolative scenarios. In al- 
ternatives to, standard percolation such as "bootstrap" 
percolationtll, one simply removes metallic sites or links 
if they do not have enough metallic neighbors. This is 
repeated until no more sites can be removed. However, 
numerical evaluation of the cluster-jnumber exponent r 
suggests no new universal behavioi£3. To generate com- 



pact configurations, Icpg-range "correlated disorder" was 
also investigated heretJ. The disorder is selected to fol- 
low the real-space correlation function C{£) — {l+i^)^'^^^, 
which falls off like i''^ for large distances £ but is well- 
behaved at £—0. Typical configurations confirm the pres- 
ence of large clusters (Fig. 2(e)). In addition, our con- 
ductivity studies for 7=0.2 (Fig. 3) are compatible with 
a neW|-universality class and i>1.31=tsp, as predicted 
beforetj. However, ta.reproduce the large exponent t~7 
found in experimentsEI anomalously small values of 7 are 
needed. The reason is that even with correlated dis- 
order the weakest links still govern the critical behav- 
ior (Fig. 2(e) shows weak links separating large clusters). 
Correlated disorder, with its statistical nature that al- 
lows for those links, is not sufficient to induce the smooth 
clusters needed to describe manganites. Naively, these 
already known results appear to rule out the possibility 
of finding an abrupt percolative transition to mimic the 
behavior of manganites. |— . 

Another alternative is the Swiss-cheese modelllZl. 
Halperin et a/O have shown that this model is equivalent 
to a random resistor network with a broad distribution 
of conductances g given by 

Pig) = il-p)S{g)+p{l-a)g-". (1) 

Here P denotes the probability that a conductance has 
values between g and g + 6g, p is the metallic fraction 
as before, and P{g) is characterized by an additional 
parameter a. The weakest bond dominates in the ef- 
fective ijesistance of the backbone. Extensive numeri- 
cal testslia support the conjecture that t=l/(l — a), for 
0.25<a<l. The previous equation becomes an equality 
for a large enough. Hence, since the data in Fig. 1(b) fol- 
lows a straight line with slope —7, it can be fitted with 
results of a simulation at a ~ 0.857. Although this nu- 
merical success is encouraging, it is still intuitively unsat- 
isfactory and does not address the first-order percolative 
transitions of Figs.l(b,c). In fact, the rapid change in p^c 
in Fig. 1(a) may also be indicative of a first-order transi- 
tion hidden by imperfections in the samples used. 

III. FIRST-ORDER PERCOLATION 

The above examples clearly show that it is highly non- 
trivial to generalize percolation to yield a first-order tran- 
sition, although it is possible to change the critical expo- 
nents. The existence of experiments (Figs.l(b,c)) with 
discontinuous transitions indicate that there must ex- 
ist yet another class of percolative transitions escaping 
our attention. We believe that the key issue missing in 
standard percolation and the related variations described 
above is a mimic of the effect of surface tension, which 
naturally arises when two phases compete. This natu- 
rally penalizes the roughness of a surface and avoids the 
weak links ubiquitous in continuous percolation. 

The simplest way to model the gregarious tendency of 
the two competing phases of manganites is to consider 
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FIG. 3: Conductivity of a 2D resistor network vs. initial 
metallic fraction, for several representative cases. N, A, C, 
and S stand for nucleated (MGE4), annealed (MGE5), cor- 
related, and standard, respectively. Both the A and S lines 
are proportional to (p — Pc)*, with t ~ 1.31. The inset shows 
finite-size scaling data for A to support universality. L is 
the linear size. Line C has a correlation power 7 = 0.2 and 
exponent t~2.1. 



the nature of the cells in the neighborhood of the cell 
under inspection. If there are enough "friendly faces" 
around, the cell will not change, otherwise it will flip to 
the competing phase. In a magnetic system, the neigh- 
borhood induces a local field that influences the center 
spin. In addition, including oxygen movement, it is im- 
possible to have an undistorted octahedron in the middle 
of a distorted JT region. The phenomenological approach 
followed hereis, based on Vichniac's majority voting cellu- 
lar automatat£l. Loosely speaking, cellular automata are 
discrete dynamical systems whose evolution is described 
by a local deterministic rule: the fate of every cell at time 
t+1 is determined by a function of the configuration of 
a few of its neighbors at time t. Cellular automata de- 
scribe an approach to equilibrium in a system dominated 
by strong, short-range interactions. Vichniact3 defines 
four different neighborhoods: Q contains only the near- 
est neighbors, V consists of Q plus the center cell, H 
consists of Q plus the second-nearest neighbors, and M 
consists of H plus the-Genter cell (see Fig. 2(c)). Unlike in 
bootstrap percolationllil, the center cell is treated on the 
same footing as its neighbors in V and M. This turns out 
to be a key property of the percolative process discussed 
here, making clusters compact and avoiding weak links. 
This is quite different from bootstrap percolation where 
the center cell is singled out. After the neighborhoods 
are defined, the update rules are as follows. Given a site 
i at a time t, let S be the number of occupied sites in 
the neighborhood TV of i. S can be any number between 
and the number of sites in the neighborhood. If S is 
greater than or equal to some threshold T, site i will be 



occupied at time t + I. Otherwise it will be epipty. The 
procedure stops when no more changes occurllj. A rule 
is specified- by a unique N and T. This is referred to 
as NGETo. For the four neighborhoods defined above, 
only 8 rules (with T close to half of the number of sites 
in N) are nontrivial, as explained below. 

The starting point for the calculations are always stan- 
dard percolation configurations, characterized by an oc- 
cupation probability p. At time t—0, p is a good mea- 
sure of the occupied fraction, within statistical fiuctua- 
tions. Under the "time evolution" provided by the cel- 
lular automaton, the occupied fraction needs no longer 
be equal to p. Rules with large (small) T tend to (in- 
crease) decrease it. For initial states with p < pc, the 
final state does not percolate. Two examples, to be dis- 
cussed in more detail in the next paragraph, are showed 
in Figs.2(f,g). Note the clear visual difference from the 
more canonical cases Figs.2(d,e). The clusters are now 
more rounded, and our intuitive goal of phenomenolog- 
ically mimicking the effect of surface tension has been 
reached. After the time evolution stops, the rejsistivity is 
measured using a numerically exact algorithmE3. 

The increase in conductivity a immediately after pc 
{— 0.5) is clearly faster in the annealed (MGE5, i.e. each 
site remains or becomes metallic if five or more of its 
nine neighbors -including itself- are metallic) case than 
in standard percolation, as Fig. 3 indicates. This is in 
agreement with our intuition that rounding effects would 
speed up percolation, mimicking better the experimental 
situation after the randomly oriented FM clusters align in 
a magnetic field. However, it is surprising that this does 
not affect the value of t, as the inset shows. Nevertheless, 
this rapid increase could be experimentally relevant. On 
the other hand, the nucleating rule (MGE4, i.e. each site 
remains or becomes metallic if four or more of its nine 
neighbors -including itself- are metallic) leads to differ- 
ent results. Its resistivity jumps from below Pc, where 
there is no connection across the system (see Fig. 2(g)), to 
the fully saturated value above Pc- To understand this, 
note that the clusters in Fig. 2(g) have very clean bound- 
aries. Adding a single metallic site immediately outside 
one of the clusters will flip all the sites touched by the 
boundary. At p close to Pc, it is likely that this apparently 
small perturbation will cause the cluster to coalesce with 
a neighboring cluster. In practice, this starts a chain re- 
action (avalanche) that only stops when the entire lattice 
turns metallic. The nucleating procedure indeed provides 
a percolative transition that is first- order. 

How common are the first-order characteristics ex- 
pected to be in the cellular automata models described 
here? Using the four (Q, V, H, M) neighborhoods defined 
above on the square lattice, there are eight non-trivial 
rules (QGE2, QGE3, VGE3, HGE4, HGE5, MGE4, 
MGE5, MGE6). By non-trivial we mean rules whose 
transition occur at a finite value of p (i.e. neither at nor 
at 1) in the infinite system limit. In all cases a careful size 
extrapolation was carried out numerically to determine 
the critical concentration. They are respectively (0.16, 
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0.83, 0.5, 0.33, 0.66, 0.25, 0.5, 0.75). Although these val- 
ues are close to the fractions (1/6, 5/6, 1/2, 1/3, 2/3, 
1/4, 1/2, 3/4), we know of no rigorous result that they 
are equal to them. In two cases not listed above (VGE2 
and VGE4) the numerical uncertainty prevented us from 
determining for sure whether the rule is trivial or not. 
To be specific, the critical p for VGE2 is very small, and 
it becomes smaller still as we increase the size of the sys- 
tem used in the simulation. We could not rule out the 
fact that p may extrapolate to zero in this case. As a 
consequence, we chose to leave those two cases out of the 
present discussion. Of the eight rules, the conductivity of 
three of them (QGE2, HGE4 and MGE4) was found to 
vary discontinuously at the transition. The conductivity 
of MGE4 is shown in Fig. 3 where the reader can con- 
vince oneself that the transition is discontinuous. QGE2 
and HGE4 look identical apart for a shift in the location 
oi Pc- Phenomenologically, the three cases behave simi- 
larly in the sense that states very close to the transition 
are unstable against flipping even just one site. Hence, 
we conclude that first-order tendencies are certainly not 
a fragile feature among the set of possible voting cellu- 
lar automata. However, far more work is needed to fully 
understand the behavior of these novel systems, with the 
present paper hopefully providing the first steps in that 
direction. The detailed analysis of the models introduced 
here may be as rewarding as the many previous studies 
of standard percolation. 



IV. DISCUSSION 

The novel discontinuous percolative transition de- 
scribed here can be applied phenomenologically J;o the 
Mn oxides. In a previous related effort, Mayr et al.til con- 
structed random resistor networks using standard per- 
colation and provided a tentative phenomenological de- 
scription of manganites. The curves pdc vs. T mimicked 
properly the results corresponding to some manganites, 
but were too smooth to accommodate the materials of 
Fig.l as special cases. The nucleating-rule procedure 
to generate metal and insulator clusters described here 
may solve this problem. The left panel of Fig. 4 con- 
tains the effective resistivity (obtained solving the Kir- 
choff equations) vs. T for configurations generated by the 
nucleating-rule, where the light (dark) regions in Fig. 2(g) 
are insulating (metallic). The individual resistivity of 
those regions is taken from experiments, and, they carry 
a T dependence, similarly as in Mayr et alH. The new 
curves are in excellent agreement with experiments (see, 
e.g.. Fig. 1(b)), including the discontinuity, and the fact 
that there is a common value of pdc at both low- and 
high-T. 

Note that the new family of percolative processes pre- 
sented here naively seems triggered only by the occupa- 
tion fraction p, not by temperature, as opposed to some of 
the results of Fig. 1. However, the ferromagnetic fraction 
in real manganites is a function of temperature: clusters 
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FIG. 4: Left panel: Resistivity pdc vs. T using the nucleating 
rule. The values of pdc fpr metalUc and insuLating bonds are 
taken from experimentsH, as incil. Followin^j, the metallic 
fraction varies linearly with temperature, p~po — a{To — T). 
We consider a—l%o, ro=100K, and label the curves by po, 
their metallic fraction at T=To. Right panel: pdc vs. T 
of the Q=f§0 Potts model (for details, see text). A cluster 
algorithmcj was used on a G4? lattice. Long simulations (sev- 
eral thousand sweeps) were required to sample both peaks 
in the probability distribution function satisfactorily. Above 
Tc there is no percolation {pdc=oo), but this can be solved 
phenomenologically by allowing for conductivity among the 
several Q clusters, as in the dashed line. 



are expected to be formed at T* (although they are not 
sufficiently large in size to create a global FM state), 
their number grows with decreasing temperature (i.e. p 
increases in the language of this paper), until the Curie 
temperature is reached. It would be incorrect to believe 
that our results are only valid for cases where the tem- 
perature is held fixed while concentrations are changed. 

Another clarification is important. Why do all man- 
ganites not have either a first- or second-order transi- 
tion? Current experimental evidence suggests that both 
cases C|an be found, although the latter tends to be fairly 
abruptu. Several aspects may be influencing these re- 
sults. First, it may occur that available single crystals are 
affected by extrinsic effects in their preparation. From 
this perspective, the current experimental results may 
vary in the future as cleaner samples are prepared. But 
there are other issues to consider. Clearly, different man- 
ganites have different degrees of intrinsic disorder, such 
as the chemical disorder emerging from having different 
ionic sizes at the A-site of the perovskite. This disorder 
will likely influence the clusters shape, and both config- 
urations Fig. 2(f) and 2(g) may emerge at different com- 
positions, leading to either rapid but continuous or first- 
order percolative processes, respectively. It is clear that 
the issue of universality of the first-order transition in 
manganites is far from clear both experimentally and the- 
oretically and more work should be devoted to its anal- 
ysis. Here, we have simply provided a plausible ratio- 
nalization of some steps involved in the CMR transition 
of the materials mentioned in Fig. 1. Whether all CMR 
materials have the same kind of first-order transitions or 
not will hopefully be clarified by future experiments. 
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Finally, although the discussion here is based on a cel- 
lular automaton, we found it reassuring that first-order 
percolative processes can also arise from well-defined 
Hamiltonians. In particular, consider the Q-states Potts 
model known to have a a first-order transition for (5>4E3. 
Simulating the model ((5=20) with MC techniques and 
using the Fortuin-Kasteleyn mappingla to construct clus- 
ters for tie various Q's, their percolative properties were 
analyzedcJ. Clusters for each Q are individually as- 
sumed metallic, and the network resistance is calculated. 
Results (Fig. 4 (right)) indeed have a pdc discontinuity 
(above at the critical temperature Tc. At T>Tc, 

no cluster percolates leading to infinite pdc, but this can 
be easily fixed assuming a finite conductivity among the 
many Q clusters. Of course, it is not at all claimed here 
that real manganites have Potts symmetry, this is just an 
illustration of concrete model Hamiltonians that presents 
the first-order transition described in the previous section 
by cellular automata rules (rather than emerging from a 
Hamiltonian approach). 



lation. The rounding of clusters induces an abrupt first- 
order transition in the model, consistent with results of 
some experiments Fig. 1(b). The present calculation sug- 
gests that many manganites do not belong to the stan- 
dard percolation class, as believed by many before, but 
correlation effects change the universality class inducing 
a faster|4iercolative process that may present hysteretic 
featuresEj. Note that the fairly "cajidom looking" config- 
uration Fig. 2(g) is self-generatedc3 by our approach, its 
phenomenological origin lies in phase competition, and 
is not pinned by explicit disorder. Thus, the first-order 
percolative transition is predicted to appear particularly 
in nominally "clean samples" and for low intrinsic chem- 
ical disorder. If explicit disorder is added to the system 
Fig. 2(g), then the transition is expected to become con- 
tinuous. More work is needed to fully clarify the fasci- 
nating physics of CMR manganites. 



CONCLUSIONS 



In summary, a phenomenological approach is proposed 
to include correlation and surface tension effects in perco- 
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